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It is nearly impossible to separate two interleaved phonebooks when held by their spines. A full 
understanding of this astonishing demonstration of solid friction in complex assemblies has remained 
elusive. In this Letter, we report on experiments with controlled booklets and show that the force 
required increases sharply with the number of sheets. A model captures the effect of the number 
of sheets, their thickness and the overlapping distance. Furthermore, the data collapse onto a self¬ 
similar master curve with one dimensionless amplification parameter. In addition to solving a long¬ 
standing familiar enigma, this model system provides a framework with which one can accurately 
measure friction forces and coefficients at low loads, and that has relevance to complex assemblies 
from the macro to the nanoscale. 


Many of us are familiar with a classical demonstration 
of the strength of friction: take two phonebooks, inter¬ 
leave their sheets and try to separate them by pulling on 
their spines. This demonstration has been carried out 
spectacularly by attempting to pull the books apart with 
people, trucks, lifting a car [Tj, and even two military 
tanks [2], only to fail and suggest that the inner friction 
between these sheets prevails. The simple explanation 
often given is that gravity provides the normal force that 
generates the tangential friction, but this hypothesis is 
easily proven to be wrong as there is no discernible dif¬ 
ference between such an experiment carried out in the 
vertical or horizontal direction. In this Letter, we study 
the force needed to separate two books as a function of 
the number of sheets, the thickness of the sheets, and 
the interleaving distance. In particular, we show that the 
force required to separate the books increases abruptly 
with the number of sheets. The strength of the system 
is due to the operator: the person, car, truck, or tank, 
amplifies any small friction arising from the normal force 
acting on the boundaries of the stack. We present a sim¬ 
ple model that captures all the data into a self-similar 
master curve. The model depends on one single dimen¬ 
sionless amplification parameter, and thus gives insight 
into the mechanisms at play in this deceivingly complex 
system. In addition to solving a long-standing familiar 
enigma related to the classical problem of friction, this 
model system provides a framework within which one can 
accurately measure friction forces and coefficients at low 
loads, and opens the way to the technologically-relevant 
engineering of friction in complex assemblies from the 
macro to the nanoscale. 

The first-known systematic studies of friction were car¬ 
ried out five centuries ago by da Vinci mm who discov- 
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ered basic rules which were later confirmed by Amon- 
tons [3j . In particular, these laws establish that the fric¬ 
tion force is independent of the contact area and pro¬ 
portional to the applied load during sliding, the propor¬ 
tionality constant being the coefficient of kinetic friction. 
Coulomb rediscovered these laws and further determined 
that, during sliding, friction is independent of the relative 
speed between the surfaces [3]. This simple set of rules, 
collectively known as Amontons-Coulomb (AC) laws, has 
been well studied in macroscopic experiments over the 
centuries. During the last decades, efforts on the micro 
and nanoscale, and towards biology [5], have resulted in 
a resurgence of activity in tribology. For example, tools 
like the surface force apparatus and the atomic force mi¬ 
croscope have fuelled experimental efforts [BHS!, as well as 
advanced theoretical treatments that go well beyond phe¬ 
nomenology pKni- Interest in the development of micro- 
electro-mechanical-systems and mechanical devices that 
operate on small length scales has driven much of this 
research. At the extreme limit, down to the nanoscale, it 
was found that the energy dissipation in friction depends 
on both electronic and phononic contributions, and that 
differences in the electron-phonon coupling between sin¬ 
gle and bilayer sheets of graphene result in variations 
in friction [12] . Furthermore, friction was probed in ex¬ 
periments on multi-walled carbon nanotubes and boron 
nitride nanotubes P3HI5]. In such investigations, the 
inner tubes could be slid out of the outer tube reveal¬ 
ing vanishingly-small molecular friction for carbon, and 
a much stronger, area-dependent, molecular friction for 
boron nitride m- Clearly, such works reveal drastic de¬ 
partures from the simple AC laws, as one approaches the 
nanoscale. Moreover, friction can be further complicated 
by a stick-slip response, as seen in many cases mm, 
including layers of paper nans. Remarkably, although 
these complex non-linear phenomena may be observed at 
micrometric scales, it is often the case on larger length 
scales that the simple AC laws with which we are so fa- 
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FIG. 1. The enigma of the two interleaved phonebooks: ex¬ 
perimental setup, and theoretical model, a) Schematic of the 
experimental setup, b) Schematic of the interleaved books in¬ 
troducing the local traction force T n exerted by the operator 
on the n th paper sheet. The two vertical lines represent the 
clamps, while the other lines represent the central lines of each 
sheet. 


miliar are valid (see reviews 0 [HI ED EMEU])- 

In this Letter, we focus on the stubborn popular 
enigma of the two interleaved phonebooks mm- In order 
to investigate these striking observations, we developed a 
well-controlled system of two identical books, both made 
up of 2 M 1 identical sheets of paper (Inacopia™) 

with width w = 12 cm, length L = 25 cm, and thickness 
e = 0.10 mm. The two books were prepared by stack¬ 
ing sheets and holding the assembly at one end with a 
rigid aluminium clamp. They were perfectly interleaved 
sheet-by-sheet, and mounted by their clamps vertically in 
a traction instrument (Adamel Lhomargy DY32). Then, 
they were separated in a vertical orientation, while the 
total traction force T was measured (see Fig.JTJa)). Since 
the experiments explored different ranges of force, three 
sensors were used with maximum forces of 10 N, 100 N, 
and 1000 N. The length of overlap between the two books 
is denoted by L — d, where d is the separation distance 
(measured with an accuracy of 10 ^tm) from the clamp 
of each book to the contact zone (see Fig. [fjb)). Consis¬ 
tent with the actual experimental parameters, we make 
the simplifying assumption that d is large compared to 
the total thickness 2 Me of one book, so that the angle 
9 n made by the n th sheet as it traverses from the clamp 
to the contact zone is small. Therefore, L — d is nearly 


identical for all sheets. Finally, the books are separated 
at constant velocity (typically 1 mm/min) and we have 
found, in accordance with the AC laws, that the velocity 
does not significantly affect the results. 

As the books are pulled apart, an initial maximum 
traction force is first reached before they start to slide 
with respect to each other. In Fig. [2j we show the raw 
total traction force T, measured during const ant-velocity 
sliding, as a function of the distance d, for seven experi¬ 
ments with M ranging from 12 to 100. As observed, the 
smaller d or the larger M, the larger the traction, and 
those dependences are highly nonlinear. Furthermore, 
the amplification of friction is far from being a small ef¬ 
fect: a single experiment spans over three decades in the 
traction force. Additionally, a tenfold increase in the 
number of sheets ( e.g. M = 12 and M = 100) induces a 
four orders of magnitude increase in the traction force. In 
the left inset of Fig. [2j we can see clear evidence of stick- 
slip mum for an experiment with M = 50. However, 
the difference between the local maxima and minima, 
which results from the difference between the coefficients 
of static and kinetic friction, is negligible in comparison 
to the global amplitude in T. We can thus neglect this 
effect, as well as the difference in the coefficients of static 
and kinetic friction. Furthermore, in the right inset of 
Fig-i we see that the experiments are reproducible and 
independent of the initial separation distance do- Fi¬ 
nally, the friction of paper can be affected by humidity 
changes [223 • AH the experiments were carried at am¬ 
bient humidity which can vary from 30 to 60 % but is 
usually closer to 45 %. The day-to-day variations are 
negligible as can be seen from the reproducibility of the 
experiments. 

These results can be explained using simple geometri¬ 
cal and mechanical arguments, for which Fig. 0b) pro¬ 
vides detailed notations. Each sheet of a given book is 
indexed by n ranging from n = 1 in the middle of the 
book to n = M at one extremity. The problem is sym¬ 
metric with respect to the central line of the book. We 
define H n = h n /d — ne/d, where h n is the shift in posi¬ 
tion of the n th sheet in the contact zone with respect to 
its position of clamping. As a consequence, the tilt angle 
8 n of the n th sheet satisfies sind n = H n /^/l + H% and 
cos 9 n = 1/ sj 1 + H%. Essentially, the tilting of each in¬ 
dividual sheet n in the intermediate region between the 
clamping and contact zones converts part of the local 
traction T n exerted on it by the operator (at point A) into 
a supplementary local normal force T n tan 6 n = H n T n 
exerted on the stack below it (at point B ). Therefore, 
according to AC laws, this leads to a self-induced addi¬ 
tional inner friction force that resists the traction: the 
more the operator pulls, the higher the frictional resis¬ 
tance. At onset of sliding, the change in traction with n 
thus reads: 

T n — T n+ 1 = 4 fxH n T n , (1) 

where /i is the coefficient of kinetic friction, and the two 
factors of 2 come from the identical contributions of the 
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FIG. 2. Measured total traction force as a function of sep¬ 
aration distance for various numbers of sheets. The number 
of sheets for each of the two identical books is 2 M, where 
M = 12,15,23,27,50,75, and 100, with M increasing as in¬ 
dicated by the dashed arrow (for M=50 and 100 only one of 
two data sets are displayed for clarity). The left, inset shows a 
zoomed region for M = 50 highlighting stick-slip friction. The 
right inset shows three experiments with M = 23, and varying 
initial separation distances do = 11,17 and, 22 mm (vertical 
arrows). 


two books and the two pages of one sheet. Finally, the 
boundary condition is given by Tm = T*, where T* is 
the unknown traction exerted on the outer sheet, and 
the total traction force is defined by T = 2 Ylk=i Tk- 
We now introduce the variable z = n/M , as well as the 
dimensionless amplification parameter a = 2peM 2 / d, 
and use a continuous description by replacing T n by T(z), 
since M 1. We thus obtain the ODE: 


T\z) + 2az T(z) =0 . (2) 

Together with the boundary condition T(l) = T*, 
we get by integration the local traction force T(z) = 
T* exp[a(l — z 2 )\. Finally, we obtain the total traction 
force T = 2 M f Q dzT(z) in the self-similar form: 


r 

2MT* 



exp(a) erf (y/a) ■ 


(3) 


This expression is remarkable for two reasons. First, it 
tends to 1 as a —i 0, which means that it represents the 
geometrical amplification gain of friction with respect to 
a simple linear collection of 2 M independent flat sheets 
with local friction T*. Secondly, it depends solely and 
almost exponentially on the amplification parameter a, 
which thus appears as the central dimensionless param¬ 
eter of this study. 

We find that the data is well described by Eq. ([3]), us¬ 
ing a single value T* = 0.01 N of the microscopic traction 
force, and with a single coefficient of kinetic friction fit 
freely to each traction experiment. We note that /i is 
not strictly constant during an experiment since the load 
varies, as explained below. Taking p to be constant is a 
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FIG. 3. Coefficient of kinetic friction as a function of the 
average effective load per page, (circles) Values of p obtained 
by fitting each experiment performed on the interleaved books 
to Eq. as a function of the average effective load (A P) de¬ 
fined in text. Since p is load-dependent, we indicate the trac¬ 
tion range over which Eq. was fit to the data by the grey 
horizontal lines, (squares) Values of p obtained independently 
using a tribometer for two sheets of the same paper. In this 
case, (M) is simply the externally-applied load, and the error 
bars correspond to the standard deviation of the measurements 
that have been repeated at least five times. The inset shows 
a fit to Eq. © (dashed line) of a traction experiment (plain 
line), for M = 50. 


simplifying assumption in the model and it captures the 
average value of p over the range in load. A typical exam¬ 
ple is shown in the inset of Fig. |3j for which an excellent 
fit is obtained over more than two decades in the traction 
force, with a constant value p = 0.73 ± 0.02. The range 
over which our simple model is applicable, with a con¬ 
stant p, is clear from the fit, and deviations are observed 
at small overlap between the books. Note that the tiny 
force T* acting on the outmost sheets is a crucial bound¬ 
ary condition. It is this finite force, corresponding to no 
more than the weight of a butterfly, that is self-amplified 
by the operator - either in a well-controlled experiment 
as we have carried out here, or when lifting a car with 
Phonebooks PQ. In our experiment, the boundary force 
T* originates from the elasticity of the paper: the outer 
sheets have a tendency to be flat and resist slightly the 
bowing induced in the contact region, thus creating a 
small normal force resulting in a small friction force. 

The coefficient of kinetic friction is a phenomenological 
quantity whose value depends on several parameters |23j . 
For metal-on-metal [24| and paper-on-paper [ 25] , p is 
known to be dependent on the normal load. Specifically, 
p increases with decreasing load and can even be higher 
than 1 for small loads, due to adhesion forces [22]. This 
fact is confirmed in Fig. [3j where we plot the coefficient 
of kinetic friction obtained from a standard tribology ex¬ 
periment between two sheets of paper, as a function of 
the applied normal load (squares). These values were 
obtained at 23 °C and with a relative humidity of 50 %. 
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FIG. 4. Self-similar master curve of the rescaled total traction 
force as a function of the amplification parameter. To produce 
this figure, each of the data sets shown in Fig. [g] (same colour 
code) was fitted to Eq. (dashed line), using T* = 0.01 N, 
and p as a free parameter (see Fig. [3]). 

Let us now compare this result with the interleaved-book 
case. As explained above, our model is based on the AC 
laws and considers a single average p for each experiment, 
that is constant irrespective of the load. This hypothe¬ 
sis is equivalent to taking a coarse-grained approach, even 
though the load increases towards the centre of the books 
and at small d. The previous assumption is sufficient to 
describe the data, as evidenced by the fit in the inset of 
Fig- § Then, we define through AC laws an effective 
load per page, N = T/(AMp), which would correspond 
to the load applied on each page if the 2M sheets were 
parallel ( i.e. all 0 n = 0). Note that, by defining the 
microscopic load J\T* = T*/{2p), Eq. ([3]) also provides 
the amplification gain N/N* in normal load. The av¬ 
erage traction force (T) of a given experiment is calcu¬ 
lated over the range of applicability of the theory, and 
used to obtain the average effective load per page (J\f). 
The result p((Af)) is shown in Fig. [3] (circles), with the 
ranges of applicability indicated by the grey horizontal 
lines. In short, a single value of p could be used to ade¬ 
quately describe each T{d) experiment with the theory, 
even though the load varied over the indicated range. 
The consistency between the simple tribology data and 
that of the interleaved books provides confidence in the 
model. Furthermore, it is noteworthy that the classic 
phonebook demonstration probes a broad range of loads, 
that typically consists of two decades, and a nearly six¬ 
fold change in the coefficient of kinetic friction. 

The self-similar Eq. ([3]) suggests that with T* = 


0.01 N, and the best fit values of p for each experiment 
(see Fig. [3]), all data should collapse when normalised ap¬ 
propriately. Indeed, in Fig. [4j we plot the rescaled total 
traction force as a function of the amplification parame¬ 
ter, and observe a single master curve for all experiments, 
with varying d and N. 

The tremendous geometrical amplification of friction 
through the interleaving of the 2 M sheets of each book 
bears a resemblance to two everyday examples. First, a 
sailor can moor a large ship simply by wrapping a rope 
around a cylinder known as a capstan [2S1127]. Note that 
its principle might be relevant to the interaction between 
DNA and a bacteriophage capsid |28j. The second ex¬ 
ample is a simple toy known as the Chinese finger trap, 
where a woven helical braid is loosely wrapped around a 
finger. The braid tightens and the finger is trapped as it 
is pulled. The trapping mechanism results from a simple 
conversion of the traction force to an orthogonal compo¬ 
nent, which enhances the load and thus the friction. This 
type of braid is applicable to sutures in surgery [21], and 
is also thought to play a role in adhesive proteins ! 30j . 

Crucial to the geometrical amplification of friction with 
interleaved books are the angles 9 n that the sheets make 
as they approach the contact region. It is through these 
angles that the traction forces result in loads perpen¬ 
dicular to the paper-paper interfaces, and thus in large 
self-created friction forces. This fact is easily verified: 
one can realise an interleaved-book system with 9 n = 0, 
by removing alternating sheets in two notepads. In such 
a case, the books can be easily pulled apart, consistent 
with the theory presented. 

Finally, it is clear from the normalised data in Fig. [4] 
that the experiments deviate from the model at very 
small loads. This is the result of assuming a coefficient 
of kinetic friction that is independent of the load. While 
it is certainly possible to include numerically an ad-hoc 
load-dependent coefficient of kinetic friction, this would 
be at the cost of simplicity and additional free param¬ 
eters. In this Letter, we have instead opted to capture 
the essential physics needed to elucidate the interleaved- 
phonebook enigma, and we reveal the key dimensionless 
parameter a = 2peM 1 2 jd of the problem. 
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